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QUANTITATIVE DETERMINATION OF OPTICAL IMPERFECTIONS 
BY MATHEMATICAL ANALYSIS OF THE FOUCAULT 
KNIFE-EDGE TEST PATTERN 

By S. Katz off 
Langley Research Center 

SUMMARY 

Linfoot developed an integral expression for the distribution of light intensity across 
the image of a lens or mirror when it is tested by the Foucault knife-edge method. This 
paper presents a fairly convenient method of inverting the linearized form of this expres- 
sion so that the errors in the figure can be determined quantitatively from measurements 
of this light-intensity distribution. Essentially, the method consists of making a Fourier 
analysis of the deviation of the light- intensity distribution from that for a perfect lens or 
mirror, and then summing the errors that correspond to the different terms of the 
Fourier series. The method is developed both for mirrors (or lenses) without central 
obscurations and for mirrors with central obscurations (Cassegrain type). 

Various examples are worked out for both types, both to show the practicability of 
the method and to show the relationships between the error in the figure and the deviations 
of the knife-edge test pattern from the normal. 

Some discussion has been included concerning application of the method to the 
testing of a large orbiting telescope, possible limitations to the basic theory of the knife- 
edge test, and the desirability of solving the same problem for the case when the errors 
are so large that the present linearized theory is either inexact or totally inapplicable. 

INTRODUCTION 

A number of laboratories, including the NASA Langley Research Center, have been 
actively involved in studying the problems associated with a large orbiting telescope. In 
order that the telescope might take full advant^e of the perfect seeing that exists outside 
the atmosphere, it should have an essentially diffraction-limited performance. Accord- 
ingly, an appreciable part of this study effort has been devoted to problems associated 
with achieving and maintainii^ such perfection. To manufacture the primary mirror for 
such a telescope, while accurately allowii^ for (or compensating for) gravity, is itself a 



formidable task. Then, even if the telescope is diffraction-limited when it goes into 
orbit, temperature nonuniformity and possibly material relaxation or creep may subse- 
quently introduce distortions. 

On the basis of such considerations, it is apparent that some method of in-orbit 
testing is desirable. If the primary mirror is made adjustable - for example, slightly 
flexible, with a system of jacks, or actuators, for adjusting the figure (see ref. 1) - it 
would be especially desirable that the test method be capable of providing quantitative 
information on the distortion that the adjustment mechanism must correct. 

The Foucault knife-edge method offers a possible approach to solving the problem 
of in-orbit testing. It is widely used in the testing of lenses and mirrors because of its 
simplicity, sensitivity, and minimal equipment requirements. The method would seem to 
be especially suited to testing the optics of a large telescope in orbit because the essen- 
tial requirement of a perfectly uniform and perfectly collimated test beam of adequate 
intensity is satisfied in outer space by the light from any sufficiently bright star. The 
usual cut-and-try method of removing any revealed imperfections, however, hardly seems 
to be applicable to a mirror in orbit, although it is possibly not completely impractical. 

In any case, a quantitative analytical approach to interpreting the knife-e(%e-test obser- 
vations seems better adapted to automation of the test and correction procedure. 

This problem — the quantitative interpretation of the appearance of the mirror 
under the knife-edge test - is the subject of the present paper. The approach starts 
from the analysis of Linfoot (ref. 2) who derived expressions for the distribution of light 
intensity in the knife-edge test pattern for the perfect mirror, and also for the imperfect 
mirror as a function of the surface distortion (more specifically, as a function of the 
intensity and phase distortion of the converging spherical wave front as it leaves the 
mirror). For a uniformly reflecting mirror with relatively slight inaccuracies in the 
figure, the corresponding change in the pattern from that for a perfect mirror is given 
by a linear integral expression involving the phase distortion of the converging wave front 
at the mirror. Solving for this phase distortion (in effect, solving for the distortion of the 
mirror surface) amounts to inverting this integral expression. 

Thus, the problem of quantitatively interpreting the knife-edge observations reduces 
to the problem of solving an integral equation. A method of solving the equation is pre- 
sented herein. In the meantime. Dr. B. E. Gatewood at Ohio State University has studied 
various other approaches under an NASA grant (ref. 3). 

The present paper reviews the method of solving the integral equation for two basic 
cases, namely (1) the mirror without a central obscuration, which will be referred to as 
a simple mirror, and (2) the mirror with a central obscuration, as in a Cassegrain 
arrar^ement. Several examples are analyzed, including the examples analyzed in 
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reference 3. Because high accuracy requires large computational effort, the accuracy in 
these examples (determined by the number of terms in the Fourier series that represents 
the solution) was limited to that svifficient to show the main characteristics of the methods. 

A few remarks have been included concerning the actual procedure of applying the 
knife-edge method to study the distortion over the entire surface of a mirror. These 
remarks are of a general nature and do not go into the important mechanical and instru- 
mentation problems that would have to be solved in order to make the knife-edge method 
practical for use with an orbiting telescope. 

SYMBOLS 

D(x,y) electric wave displacement at the (x,y) image point in the observation plane; 

normalized so that it would equal E(x,y) for perfect optics and with 
knife-edge removed 

E(x,y) electric wave displacement at (x,y) on the converging wave front just as it 
leaves the mirror; amplitude normalized to unity 

I(x) intensity (see eq. (3)) 

Iq(x) intensity for a perfect mirror 

matrices defined in the text 

local depression of mirror surface below that for a perfect mirror, in half- 
wavelengths 

vector whose elements are the values of m along a chord 
Legendre polynomial 

deviation of observed intensity from that for a perfect mirror (see eq. (4)) 

vector whose elements are the values of p along a chord 

one -half the width of gap in chord, for mirror with a central obscuration 
(For all calculated examples, R= 5/41, relative to chord or diameter 
equal to 2.) 


m(x) 



P(x) 


{p> 
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t dummy variable corresponding to x 

X,y coordinate axes 

X distance normal to knife-edge, normalized such that length of chord being 

studied is 2 

y distance parallel to knife-edge 

z distance used in defining bump on mirror surface for two examples from 

reference 3; 5x - 3 

9 dummy variable corresponding to <p 

<p = cos'^x 

Subscripts: 

a antisymmetrical 

I left side of chord (negative x) 

r right side of chord (positive x) 

s symmetrical 


OPTICAL ARRANGEMENT 

A schematic sketch of the arrangement of the optical elements is shown in sketch (a). 
The knife edge passes precisely through the focus of the mirror, and the lens, which is 
located directly behind the knife ec^e, focuses the mirror surface onto the indicated mir- 
ror image plane. This image plane might contain, for example, a raster of photosensitive 
sensors, a single line of such sensors, or possibly a single sensor mounted so that it can 
scan a line across the image. With the knife e(^e normal to the plane of the sketch, as 
shown, the scan line would be in, or parallel to, the plane of the sketch. 

Sketch (a) shows a perfectly collimated beam (as of starlight) and a parabolic mir- 
ror. In a typical laboratory test, with a spherical mirror, both the light source and the 
knife-edge are close to the center of curvature. Although such an arrangement is geomet- 
rically different from that of figure 1, both arrangements serve to produce a converging 
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Mirror- 

image 

plane 


spherical wave front which comes to a focus at the knife-edge; and the analyses of refer- 
ence 2 and of the present paper apply equally well to both cases. 

With regard to accurate placement of the knife-edge, the focus of the mirror is not 
strictly definable if the mirror is imperfect. A related fact is that, even if the mirror 
is perfect, mislocating the knife-ec^e will produce a pattern in the focal plane of the lens 
which indicates that the mirror is canted, has spherical aberration, or both (ref. 2). 

With a flexible mirror one could, presumably, then distort and rotate the mirror until it 
appears as a perfect mirror with its focus precisely on the mislocated knife-edge. How- 
ever, it would be more reasonable to try to minimize the forces exerted by the actuators 
on the mirror — that is, to set the knife-edge carefully at approximately the "best focus," 
as indicated by the general appearance of the pattern, before proceeding to analyze the 
pattern and then correct the residual distortions with the actuators. Such an approach is 
also desirable with regard to quantitative analysis of the pattern, because accurate analy- 
sis by the present linearized method does, in fact, require that the distortions be small - 
say, of the order of one-tenth of a wavelength or less. 


EQUATION FOR THE SIMPLE MIRROR 


Derivation of Linearized Equation 

General equation .- Let E(x,y) be the electric wave displacement at point (x,y) in 
the mirror surface. Let D(x,y) be the electric wave displacement at the corresponding 
image point in the observation plane (mirror-image plane, sketch (a)). In order to sim- 
plify the form of the eqviation, the measure of D(x,y) is adjusted so that if the knife-edge 
were absent along with the corresponding diffraction phenomena, D(x,y) would be the 
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same as E(x,y) for a perfect mirror. Then, with the knife-edge passing through the 
nominal focal point of the mirror, the relation between D(x,y) and E(x,y), as given by 
Linfoot, is 


2;rD(x,y) 



t - X 


( 1 ) 


where 


X distance normal to knife-edge 

y distance parallel to knife-ec^e 

t dummy variable corresponding to x 

Sketch (b) will help to clarify this equation. The circle represents the mirror (or 

its image produced by the small lens), and the knife-edge is parallel to the Y-axis. The 
form of the equation shows that 

(1) Along any chord normal to the knife edge (that is, for a fixed y), D is deter- 
mined only by E along the same chord, and thus is not affected by any other part of the 
mirror surface. This independence of each chord is a valuable simplification in ana- 
lyzing the knife-edge pattern. The parameter y will be omitted from subsequent equa- 
tions, since it will be understood that the analysis is always concerned with some partic- 
ular chord. 

(2) Since the factor dt/(t - x) in the integral is nondimensional, any convenient 
unit may be used for the linear dimensions. Any chord that is under consideration may 
thus be considered to extend from x = -1 to x = +1, as indicated in sketch (b) and in 
the limits of the integral. 
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Linearized equation .- Let E(x) = 1 for a perfect mirror. Then, if, at points x, 
the mirror surface is depressed by m(x) half-wavelei^ths, E(x) will be given by 


E(x) = e"2’rim(x) = i _ 27rim(x) - 27r^m^(x) + . . . 


When this approximation is introduced into equation (1), the equation becomes 


27tD(x) = jt 1 - 27r^m2 


1 - 27r2m2(x) + 2 f ^ dt 

t - X 

ir-2.2m(x) . f * ^ - 2.2 f ' isM 

t - X t - 


Multiplying both sides of the equation by their complex conjugates, and retaining only the 
first-order terms in m gives the following equation for the intensity: 

I(x) = 45t2|d(x)| ^ = 77^ + In^ ^^ ~ ^ ^ - 4772 m(x) In ^ ^ (3) 

The sum of the first two terms on the right side, + ln2(l - x)/(l + x), is the 
intensity Iq(x) for a perfect mirror. It is plotted in figure 1. The remaining two terms 
give the difference between the intensity for the imperfect mirror and that for a perfect 
mirror. Equation (3) is accordingly rearranged and rewritten as follows: 


I(x) - 1q(x) _ 
47T^ 


p(x) = ^ 

d_l t 


^dt - m(x) Inl^ 
X 1 + X 


The solution of this equation, where p(x) is given and m(x) is to be determined, is the 
present problem. 

Examples of p(x) Corresponding to Some 
Simple m(x) Distributions 

Substituting some simple algebraic expressions for m in the preceding equation 
will help to show some general characteristics of the relation between m(x) and p(x) . 

First, since m(x) In (1 - x)y/(l + x) is the same as ^ m(x) d^t - x), equation (4) is 


rewritten as 


p(x) = f' “W - -w 

t — X 
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Example A-1 .- If m is a constant, m(t) - m(x) is zero and hence p(x) is zero. 
In other words, a uniform phase change across the entire wave front has no discernible 
effect on the image behind the knife-edge (which fact is physically apparent, and also 
mathematically apparent since, in the derivation of the original equation, the reference 
phase was arbitrary). An alternative interpretation is that replacing the parabolic 
mirror with a smaller or larger, geometrically similar, mirror with the same focal point 
leaves the pattern unchanged when the mirror is tested with a perfectly collimated beam 
as in sketch (a) . This ambiguity is not altogether trivial, as will appear in a later sec- 
tion in which analysis of the entire mirror image is discussed. 

Example A-2 .- Let m(x) = ax; that is, let m vary linearly along the chord. Then 
p(x) = a r ^ ~ ^ dt = 2a 

d_i t - X 

Thus, for the case of a linearly varying m, p(x) is a constant. As was pointed out in 
reference 2, this case could correspond to a slight rotation of the entire mirror about a 
line through its center parallel to the knife-edge, or to a slight shift of the knife-edge, in 
the mirror focal plane, so that it no longer passes exactly through the focal point. 

Example A-3.- Let m(x) = ax2. Then 



t2-x2 
t - X 


dt = 2ax 


Other cases in which m(x) is simply an integral power of x are similarly easily 
calculated. Table I gives the first few of these and also gives the general formula. It 
will be seen from table I that if m is an even function of x, p is an odd function of x, 
and if m is an odd function of x, p is an even fimction of x. 

Table I may be inverted to derive table n, which shows m(x) corresponding to 
p(x) = x^. The coefficients in the expressions for m(x) are easily obtained by a recur- 
sive procedure, which will be given later. (See bottom of p. 12.) 


Inversion in Power Series or in Series of Polynomials 

A possible method of solving equation (4) or equation (5) is immediately apparent: 
Analyze the experimental p(x) distribution into a power series in x, 

p(x) = a + bx + cx2 + dx^ + . . . 


and the mirror distortion will be given directly by the sum of the corresponding 
terms from table II 


m(x) = -X + -x2 
2 2 





m 
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TABLE I.- p(x) FOR m(x) = x^; SIMPLE MIRROR 




A modification of this approach would be to analyze the experimental p(x) distri- 
bution into a sum of orthogonal polynomials, say the Legendre polynomials Pjj, and 
again sum the corresponding m terms. For this approach, one would need a table like 
table in, in which the p terms of table H have been combined to form the Legendre 
polynomials, and in which the corresponding m terms of table H have been combined. 


TABLE m.- m(x) FOR p(x) = P„(x); SIMPLE MIRROR 


p(x) 

m(x) 

= 1 

X 

0 

2 

X 

111 

x2 

2 

P2"i(3x2-l) 

i(i"' - 

P3 = - 2x) 


P4 = - 30x2 + 3 ) 


2\8 24 72 / 

• 

• 


A considerable amount of algebra would be involved in setting up such a table if a 
large number of Legendre polynomials is used in order to attain an acceptable accuracy, 
although the effort would hardly be prohibitive. In any case, it was finally decided that a 
similar Fourier-series approach would be somewhat more practical and generally more 
attractive. The following section describes the Fourier-series method, which is pre- 
sented as the preferred method. The material of the preceding tables will accordingly 
not be further used; it was presented because it showed most directly the relations 
between p(x) and m(x) and the nature of the proposed approach to solving the integral 
equation. 


SOLUTION IN FOURIER SERIES; SIMPLE MIRROR 

Correspondence Between p and m Expressed 
in Terms of cos”0 

In equation (5) 
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p(x) = - m(x) 

-'-I t - X 



substitute x = cos 0, t = cos 6 , where 0 is now the dummy variable corresponding 
to 0. The equation now takes the form 


p(0) = fj!!WjLjHiSi.sln9de 

•Jq cos 0 - cos 0 


( 6 ) 


Putting m(0) = 1, cos 0, cos20, cos^0, and so forth simply reproduces table I, 

X being replaced by cos 0. Inverting simply reproduces table n, similarly altered. 
These revised tables are shown here as tables IV and V. 


TABLE IV.- p(0) FOR m(0) = cos”0; SIMPLE MIRROR 


m(0) 

p(0) 

1 

0 

cos 0 

2 

COS^0 

2 cos 0 

COS^0 

2 ^cos20 + 

. 

• 


TABLE V.- m(0) FOR p(0) = cos"^0; SIMPLE MIRROR 
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Correspondence Between p and m Expressed 
in Terms of Cos n0 

One may now group the p terms of table V to form cos <p, cos 24>, by 

using the formulas (ref. 4, item 3.175) 

•N 

COS 20 = 2 cos^cp - 1 

cos 30 = 4 cos^0 - 3 cos 0 

cos 40 = 8 cos^0 - 8 cos20 + 1 

cos 50 = 16 cos®0 - 20 cos^0 + 5 cos 0 




(7) 


COS n0 = 2*^"^cos”0 - ~ 2”"^cos”"^0 

+ 2"-5cos"-^0 - "("T..4)(n - 5) 2”-’^cos”-®0 

2, 3, 

^ n(n_- 5)(n - 6)(n - 7) 2n-9^„gn-8^ _ _ 

41 


Groupii^ the correspondir^ powers of cos 0 in the m terms yields, finally, the corre- 
sponding pairs of p and m shown in table VI. The coefficients in the expressions for 
the m terms in this table are related as follows: 

1 = 1 

■ 3 ‘3 

45 5 3 3 

945 7 5 3 3 45 

428 _ 11^11^4 , 1 ^ 44 

14 175 9 7 3 5 45 3 945 

10 196 _ 1 ^ 1 ^ 1 1 ^ 4 1 s. 44 1 ^ 428 

' 467 775 11 9 3 7 45 5 945 3 14 175 
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the corresponding m(<p) is given by 


^0 ^1 ^2 / 1 \ 
m( 0 ) = cos (/> +-^ cos 20 + — f cos 30 - cos 0 j 

+ -^^cos 40 - ^ cos 20 ^ + . . . 


(It will be noted that since p(0) is defined only between 0 and n, it can be represented 
by a Fourier series containing only cosine terms.) In the example problems to be pre- 
sented herein, 12 -term series were used, for which both the programing and computing 
efforts were small in terms of modern computer technology and capability. 

The computational procedure is essentially as follows: 

(1) Locate 12 points along the mirror diameter or chord at x = cos 7.5°, 


X = cos 22.5°, X 

= cos 37.5°, 

. . ., X = cos 172.5°. The values of p at th 

(P?.5’ ^22.5’ • • 

., Pj72 5 ) 3 ^^J' 2 .i^ed in a column form the vector ^p^. 

(2) Prepare a 12 X 12 Fourier matrix [fJ as shown 



1 

1 

1 


12 

12 ' ‘ ‘ 

12 


cos 7.5° 

cos 22.5° 

cos 172.5° 


6 

6 

6 


cos 15° 

cos 45° 

cos 345° 


6 

6 

6 

[f]. 

cos 22.5° 

cos 67.5° 

cos 517.5° 


6 

6 

• • 

6 


cos 82.5° 

• • 

* • 

cos 247.5° 

cos 1897.5° 


6 

6 

6 


The product is a vector containing the 12 coefficients a^, a^, a. 2 , . . ., 

of the Fourier series. As will be seen, this step is not actually done. 
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Determining the Fourier coefficients in this manner amounts to using the trapezoid 
rule to evaluate 

2 

% = i Jq CO® 

One might suppose that a more sophisticated numerical-integration formula such as 
Simpson's one-third rule would give a more accurate result. Actually, however, in the 
case of a periodic fvinction, there is no general justification for weighting any ordinate 
more than any other ordinate. Furthermore, where a finite number of Fourier terms 
must represent a curve passing through the same number of equally spaced points (as in 
the present problem), determinii^ the coefficients by the trapezoid rule results in a curve 
that passes exactly through the given points. A proof of this interesting fact may be found 
in reference 5, section 4. 

(3) Set up the 12 X 12 matrix [m], whose elements are the 12 expressions for m 

in table VI evaluated for the 12 values of <p. Thus, the elements of the first column are 

- cos <p evaluated for 4> = 7.5°, 22.5°, . . ., 172.5°; the elements of the second column 
2 

are - cos 2<p evaluated for these same values of 4>; and so forth. 



cos 7.5° 
2 

cos 22.5° 
2 


cos 15° 
4 

cos 45° 
4 


cos 22.5° - -i cos 7.5° 

O 

4 

cos 67.5° - cos 22.5° 

O 

4 


cos 172.5° 
2 


cos 345° 
4 


cos 517.5° - cos 172.5° 
o 

4 


(4) Multiplying the vector of Fourier coefficients by the matrix [M^ finally gives 
a vector ^m^ containing the values of m at the 12 points x = cos 7.5°, x = cos 22.5°, 
. . ., X = cos 172.5°. 

Combining all these steps into one equation gives 
{m>= [m][f]^> 
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The 12 X 12 matrix [f] is thus the basic result. Multiplying any experi- 
mental p vector by |f] gives the corresponding m vector ^m^ which 
represents the mirror distortion at the 12 points. This inversion matrix (m] is 
given completely in the appendix. The number of significant figures shown is no doubt 
excessive; no effort was made to determine the optimum number of significant figures. 


Example Problems 


Three example problems were calculated by use of this inversion matrix. The 
first two were very simple ones that served mainly to verify that the basic concept of 
the method is correct and that there were no errors in the matrix. The third problem 
was the same one that was studied in reference 3. It corresponds to a narrow bump on 
the mirror, centered at x = 0.6, and of total width equal to 0.4. These three problems 
are reviewed in this section, 


Problem A-1 .- In the first problem, all the elements of the vector ^p^ were 
unity (1,1,1,. . .,1). As seen in table in or table V, the corresponding m is simply 
^ X or i cos (f). The elements of the vector ^m^ were indeed exactly ^ cos 7.5°, 

1 cos 22.5°, ^ cos 37.5°, . . ., ^ cos 172.5°. The p and m curves are shown in 
figure 2(a). 

Problem A-2.- In the second problem, p was taken as x^ or cos^(p (which is 

1 1 \ ' 
equal to — + •^ cos 2<pj. As seen in table H or table V, the corresponding m is 

or i^cos^</> - cos <p^. The calculations produced this result exactly. The 
p and m curves are shown in figure 2(b). 


The accuracies of the preceding two results simply exemplify that if p is exactly 
expressible in the form slq + a-y cos <p + ^2 cos 2<p + . . . + a^y cos 110, this procedure 
will give the corresponding values of m at the 12 points exactly. K, however, p can- 
not be satisfactonly represented with this 12-term series, the 12 X 12 inversion matrix 
will give a set of 12 m values that are correspondingly inexact. 

Problem A- 3 .- The bump used for the third problem is shown in figure 3. It is 
best described in terms of a new variable z which has its origin at the center of the 

bump (x = 0.6) and a unit that is one-fifth the unit of x. Thus, z = ^ ~ = 5x - 3. The 

0.2 

surface distortion is then given by 
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m(z) = 0 

m(z) = -(1 - z2)^ 


(z2 ^ 1) 
(z2 ^ 1) 



The minus sign is used in front of the parenthesis because, in the derivation of equa- 
tion (4), positive m was assumed to represent a depression. 

In the analytical evaluation of p for this case, some caution must be exercised in 
using equation (4) or equation (5), because the analytical expression for the bump does 
not apply outside the bump. Nevertheless, the calculation is straightforward and is not 
detailed here. The result is 


p = -2z^ + -^ z - (1 - In 

U 


1 - z 1 + X 
1 + Z 1 - X 


(z2 g 1) 


P = 


-2z3 + z - (1 - z2)^ In 

O 


1 - Z 
1 + z 


(z2 5 1) 


This p distribution is shown in figure 3. 

As before, the first step is to find the values of p at the 12 points x = cos 7.5°, 

X = cos 22.5°, . . ., X = cos 172.5°. With these values as the elements of the vector 
one can then determine the 12 elements of the vector ^m^ = The 12 

m values so obtained are plotted in figure 4. A curve faired through these points 
would show the bump, but not very accurately; in particular, the bump would be too low, 
although it would be in about the right location. The poor accuracy reflects the fact that 
only three points fall in the region of the main irregularity of the p curve, which are 
inadequate to define the p curve in the region (the region of the bump). It is apparent 
that if the mirror distortions can be as sharp as this bump, more than 12 points would 
have to be used (with correspondingly larger [m], [f], and [M][?i] matrices) if satis- 
factory accuracy is to be achieved. In the methods that were studied in reference 2, 

40 points were used, with excellent results. 

Although this bump probably represents an imrealistically sharp and narrow irregu- 
larity, it was decided to pursue the matter further by doubling the number of points. A 
24 X 24 |jMi] [f]] matrix was derived by extending the procedure used for deriving the 
12 X 12 matrix; and a ^p^ vector was determined by finding the values of p at 
X = cos 3.75°, X = cos 11.25°, x = cos 18.75°, . . ., x = cos 176.25°. The resulting 
m values, shown in the lower part of figure 4, are in very good agreement with the true 
m curve. The 24 X 24 [m][^ matrix is presented in the appendix. 

The fact that the entire computed m curve falls below the original m curve 
reflects the fact that the average value of a sum of cos n<p terms (averaged with 
respect to 0, in the range between 0 and ■n) is zero. As explained previously, such a 
downward shift of the entire m curve by a constant is of no physical significance. 
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EQUATION FOR MIRROR WITH CENTRAL OBSCURATION 
Derivation of Linearized Equation 

For a mirror with a central obscuration (Cassegrain t 3 q)e), one may again, for 
simplification of the arithmetic, consider the two ends of the diameter or of a chord to be 
at -1 and +1, and the breaks to be at -R and +R, as indicated in sketch (c). With 



Sketch (c) 


E in equation (1) equal to zero in the region between -R and +R, equation (2) now 
takes the following form: 



Multiplying both sides of the equation by their complex conjugates and retaining 
only the first-order terms in m gives the following equation for the intensity: 

I(x) = 477^ |d(x)| ^ 



2^jLiX _4^2j^(x)in 
X - R 1 + X 


X + R 1 - X 
X - R 1 + X 


( 8 ) 
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Let Io(x) represent the intensity distribution for a perfect mirror (ref. 2), 


Ti^ + In'^ 


X + R 1 - X 
X - R 1 + X 


Then 


I(x) - Iq(x) 

477 ^ 



m(t) 

t - X 


dt + 


m(t) 
Jr t - X 


dt - m(x) In 


X + R 1 - X 
X - R 1 + X 


( 9 ) 


The problem is to determine m(x) along the two parts of the chord (that is, between -1 
and -R, and between R and 1) when p(x) is given along the same two parts. In the 
following section, some examples will be shown of p(x) for several simple functions 
m(x); and then the general method of solution will be developed, as before, in cosine 
series. 


Examples of p(x) Corresponding to Some 
Simple m(x) Distributions 

Example B-1 .- If m(x) is a constant, equation (9) gives p(x) = 0, as it should. 
Example B-2 .- If m(x) = p where a is a constant. 


P(x) 


a + a f 

dt _ ^ In 

X + R 1 - X 

J_1 t(t - x) J 

R t(t - x) X 

X - R 1 + X 


= 0 


Thus, m(x) = Constant and m(x) = ^ are both solutions of the homogeneous 


Constant 


is a solution was first 


integral equation, p(x) = 0. ^The fact that m(x) 

recognized by Dr. Gatewood.^ The ambiguity introduced by these two solutions, and 
methods of overcoming it, are discussed in a later section. 


Example B-3.- If m(x) = x. 


p(x) = C + C - X In 

t - X Jr t - X 


X + R 1 - X 


X - R 1 + X 


= 2(1 - R) 


If m(x) is x^, x^, . . ., then p(x) is similarly easily evaluated. The first few 

p terms for these cases are listed in table Vn, and four are plotted, together with the 

corresponding m terms, in figure 4. For all calculated examples, R was taken 
as 5/41, the same value that was used in reference 3. It will be seen that, in general, 

if m(x) is symmetrical, p(x) is antisymmetrical, and if m(x) is antisymmetrical, 

p is symmetrical. 
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TABLE Vn.- p(x) FOR m(x) = x*^; MIRROR 
WITH CENTRAL OBSCURATION 



For the following three examples, and for the remainder of the discussion of the 
Cassegrain- mirror problem, it will be convenient to differentiate between the left- and 
right-hand sides of the mirror by means of the subscripts I and r. 

Example B-4 .- If mj,(x) = +1 and m^(x) = -1, 


P^(x) = 
Pr(x) = 


_ dt ^ 

'^-1 t - X Jr t - : 


+ In 


X + R 1 - X 
X - R 1 + X 


= 2 In 


1 - X 
X - R 


dt ^ r 

^ - In 

X + R 1 - X 

= -2 In 

X + R 

J_1 t - X o; 

R t - X 

X - R 1 + X 

1 + X 


It will be seen that these two expressions have the same value at symmetrically located 
points ^that is, Pj,(x) = p^(-x)), as they should, since m is antisymmetrical. 

Example B-5 .- If mj.(x) = x and m^(x) = -x. 



t dt 

t - X 



t dt 

t - X 


+ X In 


X + R 1 - X 
X - R 1 + X 


= 2x In 

X - R 


Similarly, 


Pr(x) 


-2x In 


X + R 
x+ 1 
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In this case, Pr(x) = -p^(-x); that is, p is anti symmetrical, corresponding to the fact 
that m is symmetrical. 

Example B-6 .- An effort was made to find a function m(x) corresponding to 
p^ = -1, Pj. = +1. No simple solution to this problem was found. However, by trial, a 
symmetrical fimction 


mj,(x) = m^(-x) = 0.94873|x| - 0.96455\/|x| - 


0.51652 



( 10 ) 


was found that when substituted into equation (9), produced the desired p function to an 
accuracy of about 1 percent. Equation (10), of course, applies only for the case R = 5/41. 


SOLUTION IN FOURIER SERIES; MIRROR WITH 
CENTRAL OBSCURATION 


Decomposition of p(x) and m(x) Into Symmetrical 
and Antisymmetrical Parts 

A function such as m(x) extending between -1 and +1 can always be decomposed 
into a symmetrical and an antisymmetrical part, where the symmetrical part is 


ms(x) = 


mj.(x) + m^(-x) 


and the antisymmetrical part is 


m. 


^ ^ _ mr(x) - m^(-x) 
- 


where the subscripts s and a indicate "symmetrical" and "antisymmetrical," respec- 
tively. If the corresponding p(x) is similarly decomposed into symmetrical and anti- 
symmetrical parts, Pg and Pg^, inverting Pg will give m^, and inverting Pg^ will 
give mg; and the total solution will be, by the two preceding equations, 

m_(x) = m_(x) + m„(x) 


m^(x) = mg(x) - ma(x) 
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Solution for Antisymmetrical p, Symmetrical m 

Development in power series in cos <p .- Given that mp(x) = m^(-x), the first inte- 
gral in equation (9) can be rewritten by substituting t = -t, as follows: 


r*-R m;(t) pR m,(-t) 

\ — i— dt = \ 

t - X dj t + X 


dt = 


•Jr t + X 


Then equation (9) becomes 


PrW = ”rW(j4- 


t + X/ 


dt - mj.(x) In 


X + R 1 - X 
X - R 1 + X 


1 mj.(t) - mr(x) 


= 2x 


r 

^R 


t2-x2 


( 11 ) 


With the substitution x = cos <p, t = cos 9, the results shown in table V TTT are readily 
derived for mj.(0) = cos^cf), cos^4}, cos^<p, .... The subscript r is not used for the 
column headings in the table, since the table actually applies to both the left and right 
sides; however, the subscripts s and a are used, since they are pertinent. 


TABLE Vni.- p(0) FOR m(0) = cos2«0; MIRROR 
WITH CENTRAL OBSCURATION 


mg(0) 

Pa(^) 

COS^(p 



2(1 - R) cos 0 

COS^0 


2 

~~ g “ 

(1 - R) COS^0 + ^ ~ ^ cos 0 

_ O _ 

COS®0 

2 

(1 - R) cos^cf) + ^ ^ cos^cp + 1 " 

_ 3 5 


cos <P 


This table may be inverted just as table IV was inverted to form table V. However, 
because of the various powers of R that occur in the coefficients, the procedure that 
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was used for the simple mirror would have involved a burdensome amoimt of algebra. 
Instead, the table was set up in a recursive form, in which each line is determined from 
the preceding lines. This formulation is shown in table IX. 

By use of the formulas for cos n.9 (cos 30 = 4 cos^0 - 3 cos 9, cos 59 = 16 cos®0 
- 20 cos^0 + 5 cos 9, and so forth; see equation (7) for the general formula), the lines of 
table IX can be grouped to give table X. If now the antisymmetrical part of p(4>) is 
analyzed into a Fourier series of the form a.^ cos 4> + cos 30 + ag cos 50 + . . ., the 
correspondii^ symmetrical part of m(0) is found directly by using table X. 

TABLE IX.- m(0) FOR p(0) = cos2«+l0; MIRROR 
WITH CENTRAL OBSCURATION 


Pa^^) 

cos 0 

COS^0 

COS^0 


mg(0) 

^ cos^0 = 


2(1 


2(1 - R) 

■^|^os40 - |(1 - r 3)M^J = Mj 


2(1 


cos 60 - 1(1 - r3)M« - 1(1 - r 5)M, = Mp 

- R) |_ 3 ^5 _j ‘ 


TABLE X.- m(0) FOR p(0) = cos(2n +1)0; MIRROR 
WITH CENTRAL OBSCURATION 


Pa(^) 

mg(0) 

cos 0 
cos 30 
cos 50 

Ml 

4Mg - 3Mj 
16Mg - 2 OM 3 + 5Mj 
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The procedure is essentially similar to what was described for the simple mirror, 
except that the m terms remain as powers of cos 4>. It is described in more detail 
after the following remark. 

Remark on the Fourier a nalysis.- The function p^^ is defined only between x = R 
and X = 1 (or between <p = cos“1r and 0 = 0). In order to determine the coefficients 
of its Fourier series in the usual manner, however, Pg^ will also have to be defined for 
values of <f> that fall within the obscuration — that is, for <p between w/2 and 
cos~1r. If a large number of terms are used in the Fourier series, this additional sec- 
tion of the Pg distribution may, presumably, be chosen arbitrarily, since the series is 
required to represent Pg only on the actual mirror. However, it would be reasonable 
to choose this additional section so that it forms a smooth continuation of the given Pg 
curve, in order that the Pg curve be represented as accurately as possible by a given 
number of Fourier terms. 

In the problems to be discussed later, nine-term Fourier series were used to 
represent Pg and Pg, and of the nine points (at 0 = 5°, 15°, 25°, . . ., 85°) that were 
used in the Fourier analysis, only the innermost one, at 0 = 85°, was within the central 

obscuration ^between 0 = ^ 0 = cos'^R^ For all but one of the problems, Pg 

and Pg were given by continuous analytical expressions; hence, the values of these 
expressions at 0 = 85° were used for the innermost point as the obvious best choice, 
with excellent results. In the final problem. Pa = 1 (that is, Pj. = 1, p^ = -1), the 
p curve has too violent a break between the right and left side to yield a very smooth 
m curve. This problem is discussed separately, and a simple method of handling such 
cases is suggested. 

Inversion matrix for a ntisymmetrical p, symmetrical m.- As in the procedure 
for the simple mirror, the first step is to determine the coefficients for a cosine series 
that represents Pg. A nine-term series was used in the present work. The procedure 
is essentially as follows: 

(1) Locate nine points along the mirror semidiameter or semichord at x = cos 5°, 

X = cos 15°, X = cos 25°, . . ., x = cos 85°. The values of Pg at these nine points, 
arranged in a column, is the vector ^Pg^. 

(2) Prepare a 9 X 9 Fourier-analysis matrix |f^ . The product F J{Pa^ is a 
vector whose elements are the nine desired Fourier coefficients. As with the simple 
mirror, this step was not actually performed. 
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III Ml 


III! ■■Ill nil II I II 


III mill II II II I I III 



pa]- 


Mo = 


N- 


cos 5° 

cos 15° 

cos 25° 

cos 85° 

4.5 

4.5 

♦ • • 

4.5 

4.5 

cos 15° 

cos 45° 

cos 75° 

cos 255° 

4.5 

4.5 

4.5 * ‘ ' 

4.5 

cos 25° 

cos 75° 

cos 125° 

cos 425° 

4.5 

• 

4.5 

• • • 

4.5 

• • 

■4.5 

• 

• 

cos 85° 

cos 255° 

• • 

• • 

cos 425° . . . 

cos 1445° 

4.5 

4.5 

4.5 

4.5 

se of the formxilas in table DC, calculate M-. = — ; — — - 

1 2(1 - R) 

■jcos^0 ~ etc., for the nine values of 4> 

nge in a 9 X 9 matrix [Mg^ as shown. 

cos^4>, 

and the given 

Mj(5°) 

Mg(5°) 

M5(5°) 

Mi7(5°) 

Mj(15°) 

M3(15°) 

M5(15°) 

Mj,7(15°) 

Mj(25°) 

M3(25°) 

Mg(25°) . . . 

• • 

Mj7(25°) 

Mj(85°) 

M3(85°) 

• • 

• • 

Mg(85°) 

M^7(85°) 


(4) Set up a 9 X 9 matrix jCg^ whose elements are the coefficients in table X. The 
coefficients in each line of the table, in reversed order, form the elements of one of the 
columns in jc This matrix is given completely. 


25 




1 

-3 

5 

-7 

9 

-11 

13 

-15 

17 


0 

4 

-20 

56 

-120 

220 

-364 

560 

-816 


0 

0 

16 

-112 

432 

-1232 

2912 

-6048 

11 424 


0 

0 

0 

64 

-576 

2816 

-9984 

28 800 

-71 808 

[Cs] = 

0 

0 

0 

0 

256 

-2816 

16 640 

-70 400 

239 360 


0 

0 

0 

0 

0 

1024 

-13 312 

92 160 

-452 608 


0 

0 

0 

0 

0 

0 

4096 

-61 440 

487 424 


0 

0 

0 

0 

0 

0 

0 

16 384 

-278 528 


0 

0 

0 

0 

0 

0 

0 

0 

65 536 


It will be seen that the product [M^ [c^ is a matrix whose ij element is the 
contribution of the (unit) jth cosine at the ith point. For example, the element at the 
third row and fifth column is the value of mg at 0 = 25° (the third point) corre- 
sponding to Pg^ = cos 94 > (the fifth cosine, see table X), Multiplying the vector of 
Fourier coefficients, given by |F^<^Pg^, by this matrix jc J will now give the 

nine elements of (that is, the values of mg at 0 = 5°, 15°, 25°, . . ., 85°). 

Thus, 


{“s} [P ®] a]{Pa} 

The complete inversion matrix is given in the appendix. Multiplying the 

vector by this matrix will give {™g^- 


Solution for Symmetrical p, Antisymmetrical m 
Given that mp(x) = -m^(-x), the same procedure as before provides the equation 

Pr W ■ "■rW - “rW 1” 




mr(t) - mr(x) 

* . + 


R L t - X 

1 tnij,(t) - xmp(x) 


t + X 


•^R 


t2 - x2 


dt 


( 12 ) 
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Constant 


is a solution of the 


^Incidentally, it is obvious from this equation that m(x) = 

homogeneous equation p(x) = 0. (See example B-2.)) Letting mj.(x) be x, x^, 
x^, . . . (or letting mj.(0) be cos <t>, cos^<f>, cos^0, . . .) and proceeding as before 
yields eventually the relations shown in table XI, which corresponds to table DC, pre- 
sented previously for the case of antisymmetrical p, symmetrical m. 


TABLE XI.- m(0) FOR p(0) = cos^H; MIRROR 
WITH CENTRAL OBSCURATION 


1 

COS^(p 

cos'^0 



There is an obvious relationship between the M terms in this table and those in 
table DC; namely, Mj = cos 0, Mg = Mg cos 0, Mg = M^ cos <p, ... . 

The procedure now follows the same pattern as before. The steps are 

(1) Find the values of the symmetrical part of p, that is, Pg, at the nine points 
0 = 5®, 15°, 25°, . . ., 85°, and arrange them in a column vector 

(2) Prepare a 9 X 9 Fourier matrix Ps] as shown, such that, when ^Pg^ is 
multiplied by it, the components of the product vector will be the coefficients of the series 
Uq + ag cos 20 + a^ cos 40 + . . . + a^^g cos 160 that equals Pg at these nine points. 
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1 

9 

cos 10° 

1 

9 

cos 30° 


4.5 

cos 20° 

4.5 

cos 60° 


4.5 

4.5 


cos 80° 

cos 240° 

- 

4.5 

4.5 

(3) Set up a 9 X 9 matrix |m^ as 


Mo(5°) 

M2(5°) 


Mo(15°) 

M2(15°) 

[“a] = 

Mq(25°) 

M2(25°) 


Mq(85°) 

M2(85°) 


1 


1 

9 

• • • 

9 

cos 50° 


cos 170° 

4.5 

• • • 

4.5 

cos 100° 


cos 340° 

4.5 

• • ♦ 

4.5 


cos 400° 
4.5 

. . . 

cos 1360° 
4.5 

shown, by using the formulas of table XI 

M4(5°) 

. . . 

Mig(5°) 

M4(15°) 

. . . 

Mi6(15°) 

M4(25°) 

• • • 

Mi6(25°) 

M4(85°) 

• 

Mjg(85°) 


As mentioned previously, every term in the matrix [MgJ is cos 0 times the corre- 
spending term in this matrix |mJ. 

(4) Set up a 9 X 9 matrix jC J whose elements are the coefficients in the expan- 
sion of cos n0 in powers of cos 0 (cos 0=1, cos 2<j> = 2 cos^0 - 1, 
cos 40 = 8 cos^0 - 8 cos^0 + 1, and so forth. (See eq. (7).) This matrix is here 
given completely. 
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1 

-1 

1 

-1 

1 

-1 

1 

-1 

1 


0 

2 

-8 

18 

-32 

50 

-72 

98 

-128 


0 

0 

8 

-48 

160 

,-400 

840 

-1568 

2688 


0 

0 

0 

32 

-256 

1120 

-3584 

9408 

-21 504 

[ Ca ] = 

0 

0 

0 

0 

128 

-1280 

6912 

-26 880 

84 480 


0 

0 

0 

0 

0 

512 

-6144 

39 424 

-180 224 


0 

0 

0 

0 

0 

0 

2048 

-28 672 

212 992 


0 

0 

0 

0 

0 

0 

0 

8192 

-131 072 


0 

0 

0 

0 

0 

0 

0 

0 

32 768 


Finally, 

{“a}" [pa] jjsJ{Ps} 

The complete inversion matrix [m J ^ J [f J is given in the appendix. Multiplying 
vector {Pg^ by this matrix will give {^a^* 

Example Problems 

Problems B-1 to B-4 .- Four simple cases were calculated in order to verify the 
general correctness of the approach and also to verify the numerical accuracy of the two 
inversion matrices The assumed p functions, which 

were either symmetrical or anti symmetrical, are listed in the following table, together 
with the corresponding m functions: 


P(x) 

m(x) 

1 

X 


2(1 - R) 

2(1 - R)x 

x2 

1^1 - R)x2 + ^"^3^ 

x3 

1^1 - R)x3 + 

x4 
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In these four cases, when x is replaced by cos <p, both m(0) and p{<p) are express- 
ible in terms of a constant, cos 4>, cos 20, cos 30, and cos 40. Hence, they are 
expressible exactly by the nine-term Fourier series, and the calculated values of m 
should be exactly correct, as in fact they were. These four examples are shown in the 
four upper parts of figure 4. In all calculations, R was taken as 5/41, the same value 
that was used in reference 3. 

Problem B-5 .- In this problem, for antisymmetrical p, Pj. was assumed to be 

1 -I- 

2x In ^ . This Pj. corresponds to the case m^. = x, m^ = -x, which was discussed 

earlier. Unlike the p terms of the preceding four cases, this Pj. is not expressible 
exactly in terms of a few cos n0 terms. However, its plot is a very smooth curve, 
which approaches the origin smoothly (fig. 4(e)). The computed values of m^. fell very 
accurately on a straight line (the line designated ”m" in figure 4(e)). Accordingly, the 
points are not shown. The ticks on the curves of figure 4(e) (which have been omitted 
from the curves of figures 4(a) to 4(d) because of the small scale) indicate the x loca- 
tions of the p and m points. 

It may be noted that the calculated mj. line is somewhat below the line mj. = x. 

As noted previously, m(x) can be determined (and need be defined) only to within a con- 
stant. Accordingly, the difference is of no consequence. 

Problem B-6 .- In this problem, also for antisymmetrical p, Pj. was assumed to 
be a constant, 1 (with p^ = -1). As mentioned earlier, a reasonably accurate, but not 
exact, solution to this problem has been foimd in the form of the following algebraic 
expression: 


mp(x) = m^(-x) = 0.94873|x| - 0.96455^/|x| - 


0 . 51652 - 4 = 

i/N 


( 10 ) 


The nine values of mj. foimd by the Fourier analysis method are shown by the circled 
points in figure 5, The upper line is a plot of the preceding algebraic expression for m^, 
with a constant added in order to provide a best fit. The scatter of the points about the 
line reflects the inadequacy of the nine -term cosine series in representing a Pj. curve 
that has a fairly abrupt irregularity. Nevertheless, on the whole, the fit seems to be 
fairly good. 


A second calculation was made in which the sharpness of the irregularity was some- 
what alleviated by slightly reducing the value of pj. for the one point that falls within the 
obscuration. For this case, the elements of the vector ^p^ were taken as 1,1, 1,1, 1,1,1, 
1, and 0.9 (see the sketch in fig. 5). The calculated result is shown by the square symbols 
in the figure; and the line through them is a plot of the algebraic expression for mj- with 
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a different constant added to provide a best fit. On the whole, the fit seems to be some- 
what improved, with reduced scatter of the points about the line. 

A third calculation was made in which the elements of were taken as 1, 1, 1, 1, 

1,1, 1,1,1, and 0.7. (See the sketch in fig. 5.) Although the sketch now shows a definite 
deviation of the p curve from the straight line in the region just outside the obscuration, 
the calculated mj. points (triangles) showed very good agreement with the true curve 
(the bottom curve of fig. 5). 

It appears from this example that the present method will provide a reasonably 
satisfactory solution even in the case of a fairly abrupt change in the value of p across 
the obscuration. Judiciously choosing the p values within the obscuration can reduce 
the sharpness of the break and improve the results. 

Increasii^ the number of points (and the number of cosine terms in the series) will 
improve the accuracy. Another approach is to decompose Pg^ into two parts, one of 
which is smooth and the other of which is like this problem, and use equation (10) as the 
solution of the latter part. 

Problem B-7 .- The final problem is the same as that used in reference 3. It is a 
bump that is identical in size, shape, and radial location with the bump of problem A- 3 
for the simple mirror. It is shown together with the corresponding p function in fig- 
ure 6. The p function in this case is given by the equations 


p(x) = -2z^ + -^z - (1 - z2) In 

O 


1-zl+xx-R 


1 + z 1-xx + R 


p(x) = -2z^ + ^z - (1 - z2)^ In 


z - 1 
z + 1 


(z2 s 1) 
(z2 g 1) 


where 


z = 2.^ = 5x - 3 

0.2 

The first step in the analysis is to decompose p into a symmetrical and an anti- 
symmetrical part. The method has been described previously; table XU shows the 
arithmetic. The table also shows the nine values of mg calculated from the nine values 
of Pg and the nine values of mg calculated from the nine values of Pg. The last two 
columns list the m values that represent the final solution. 

This solution is plotted in figure 6. It is generally closer to the true solution (the 
original bump) than was the first result for the simple mirror, where the calculation used 
only six points on each side; but it is not quite as close as was the second result for the 
simple mirror, where 12 points were used on each side. 
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TABLE xn.- STEPS IN DECOMPOSING p AND COMPOSING m; EXAMPLE B-7 


<t>, deg 

X 

(cos <l>) 

Pr(x) 

P^(-x) 


Pa 

“a 

“s 

iflj.(x) 

(mg + m^) 

m^(-x) 
(“s - “a) 

5 

0.99619 

0.55995 

-0.13395 

0.21300 

0.34695 

-0.00192 

-0.00632 

-0.00824 

-0.00440 

15 

.96593 

.61076 

-.13655 

.23710 

.37366 

.00820 

.00728 

.01548 

-.00092 

25 

.90631 

.74645 

-.14199 

.30223 

.44422 

-.01098 

-.01781 

-.02879 

-.00683 

35 

.81915 

1.14548 

-.15075 

.49737 

.64812 

.01144 

.01226 

.02370 

.00082 

45 

.70711 

1,36675 

-.16376 

.60149 

.76525 

-.24742 

-.25694 

-.50436 

-.00952 

55 

.57358 

-1.53599 

-.18254 

-.85926 

-.67672 

-.45945 

-.45735 

-.91680 

.00210 

65 

.42262 

-1.70324 

-.20977 

-.95651 

-.74679 

-.06215 

-.07554 

-.13769 

-.01339 

75 

.25882 

-.66017 

-.25036 

-.45527 

-.20490 

.00653 

.00643 

.01296 

-.00010 

85 

.08716 

-.42551 

-.31433 

-.36992 

-.05559 

.02349 

-.00122 

.02227 

-.02471 


The additive constant .- In the earlier case of the simple mirror, the solutions 
for m were always obtained in such a form that the average value of m (averaged 
with respect to 0) was zero. That is, whatever may have been the original m function, 
when it was rederived from the corresponding p function, an additive constant had been 
introduced such that its average value was zero. 

In the present solution for the Cassegrain-type mirror, a different condition is 
imposed, namely, that the solution pass through the origin. That is, whatever may have 
been the original m fimction, when it is rederived from the corresponding p function, 
an additive constant has been introduced such that m = 0 at x = 0. The reason will be 
apparent on referring to table IX or table XT, where it can be seen that all the M terms 
contain cos 0 as a factor and are therefore zero at 0 = 90° (x = 0). In further clari- 
fication of this matter, consider, for example, Mj in table IX, which is listed as 

?: — cos^0. Since a constant may be added without affecting the relation between p 

1 f 1\ cos 20 

and m, Mi might also have been given as — — = cos 0 - 4 = , ^ , in which 

’ 1 ^ ^ 2(1-R)\ 2/ 4(1 - R)’ 

case its average value would have been zero. In the case of the simple mirror, all the 

M terms were actually written in this manner - that is, in terms of cos n0 - so that 

the average value of every M was zero. In the case of the Cassegrain-type mirror, 

however, the M terms are all in terms of cos*^0, and the average is zero only when n 

is odd (as in table XI) . 


NOTES ON APPLICATION OF METHOD 


Survey of Complete Mirror Surface 

Simple mirror .- The analysis in this paper has been concerned with the mirror- 
surface distortion aloi^ the diameter or any chord that is normal to the knife-edge. In 
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order to survey the entire mirror surface, many chords and/ or diameters must be 
examined. Some care must be exercised, however, in selectii^ the set of chords and 
diameters. 

Consider first the case of the simple mirror. With a fixed knife-e<%e setting, one 
could study all chords normal to the knife-edge and thereby cover the entire mirror. 
However, even if all these chords appeared to be perfect, one could not conclude that the 
mirror was perfect. As previously pointed out, one could conclude only that all these 
chords were perfect parabolas with their foci on the knife-ecfee. Additional data would be 
required in order to establish that the mirror is a perfect paraboloid of revolution. For 
example, one could rotate the knife edge throv^h 90° and examine the diameter, which 
now intersects the complete set of chords that were tested with the previous knife-e(^e 
setting. (See sketch (d) .) If this diameter also turns out to be a parabola with its focus 
on the knife edge, the perfection of the mirror would be established. Or, if the mirror 
has imperfections, the data for the set of parallel chords combined with the data for this 
diameter are sufficient to analyze the mirror svirface completely. 




An alternative approach is to rotate the knife-edge through a series of small steps 
and study only the diameter for each knife-edge setting. (See sketch (e).) If all the 
diameters are foimd to be perfect parabolas with the same focal point, the fact that all 
the vertexes coincide (since all the parabolas pass through the center point of the mirror) 
would establish that the mirror is a perfect paraboloid of revolution. Or, if any one of 
the diameters shows an imperfection, the imperfection may be evaluated quantitatively 
relative to an exact parabola that passes through the common vertex (the center point of 
the mirror), which serves as the common reference point. This method seems prefer- 
able to that of the preceding paragraph, because it can use a single row of sensors with 
fixed spacing. 
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Cassegrain mirror.- For the Cassegrain mirror, the fact that the central area of 
the mirror is obscured makes it impossible to use the center point of the mirror as a 
common reference point. Accordingly, the method of the preceding paragraph is not 
sufficient to show whether the mirror is a perfect paraboloid of revolution - every 
diameter could be a perfect parabola with its focus on the knife e(%e, yet the vertexes of 
all these parabolas need not coincide. 

Some point located outside the obscuration will have to serve as a reference point, 
and all parabolas either (1) must pass through this point or (2) must be more or less 
indirectly referable to this point. The left-hand sketch (sketch (f)) shows the first 
approach, in which all the parabolas pass through the common reference point, P. It will 
be recalled, however, that any analysis of a parabola that passes through the obscuration 

may be in error by Am(x) = where the constant is arbitrary. Accordingly, 

this approach seems to be not quite adequate, since some of the parabolas pass through 
the obscuration, as shown in sketch (f). 



The method indicated in the middle sketch (sketch (g)), in which the mirror is 
covered with a set of equal unbroken parabolas, avoids this difficulty. Every one of the 
indicated family of parabolas cuts at least one of the two parabolas AB and CD through P, 
and can thus be referred to P. This procedure has the additional advantage that all the 
test chords have the same length, so that a sir^le row of sensors with fixed spacing can 
be used for the entire survey. 


The method of the preceding paragraph avoids having to deal with the case of a 
broken parabola - that is, a parabola that passes through the obscuration. The possible 

error Am(x) = , however, may generally exist to only a negligible extent; and 

it may be unrealistic to forego using any broken parabolas for a reason that could be more 
theoretical than real. In any case, a combination of broken and unbroken parabolas can 
be arranged for which the analysis of the mirror surface would provide a unique answer. 
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Such a combination is indicated in the right-hand sketch (sketch (h)). For the three 
unbroken parabolas the errors, if any, can be determined at every point relative to the 
point P, as discussed in the preceding paragraph. Specifically, the errors at points Q 
and R are known, and hence the difference between the errors at points Q and R is known. 
If analysis of the broken parabola AB provides a different value for this difference, the 

discrepancy determines the constant in the extraneous solution Am(x) = that 

has to be removed. The entire mirror surface can thus be analyzed by means of broken 
parabolas passing through the center, all checked against the initial set of three unbroken 
parabolas. 


Optical Theory 

Cassegrain arrangement .- A mirror with a central obscuration would typically be 
used with a secondary mirror near its focus; the secondary mirror reflects the con- 
verging beam back through the hole to the instrument chambers. The effect of the sec- 
ondary mirror has not been considered in this study, where it has been assumed that the 
knife-edge is simply at the focus of the primary mirror. Some study is needed to assess 
the effect of the secondary mirror on the applicability of the present theory to deter- 
mining distortions of the primary mirror. Study is also needed to determine whether the 
effects of small errors in the secondary mirror can be analyzed as if they were due to 
errors in the primary mirror, and the errors then compensated for by distorting the 
primary mirror. (From a cursory examination of this question, it would appear that, in 
general, complete compensation is unlikely. Although the compensation may cause the 
light from an on-axis star to come to a sharp focus, the compensation would hardly apply 
for the light from off-axis stars; hence, in general, the width of the field would be reduced, 
unless the error in the secondary mirror was very small and gradual.) 

Diffraction theory .- There are some approximations in the development in refer- 
ence 1; however, they seem to be of quite minor significance if the total ratio is fairly 
large. One assumption, of uncertain significance, is that the knife-edge does not affect 
the part of the wave front that does not impinge on it. This assumption has been of con- 
cern in the past where it has been applied in diffraction theory; and its validity may be 
especially questionable in the present case, where the entire Airy disk is only a few wave- 
lengths in diameter. In any case, some experimental quantitative verification of the basic 
theory of reference 2 is desirable. 

Decomposition of p and m for the Simple Mirror 

In the analysis for the mirror with the central obscuration, p was decomposed 
into a symmetric part Pg and an antisymmetric part Pg^, and m was then foimd in the 
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corresponding two parts which then had to be combined. This procedure had the advan- 
tage that two 9x9 matrices had to be developed and used instead of one 18 x 18 matrix, 
which would have been more troublesome. It should be pointed out that this procedure is 
directly adaptable to the simple mirror case. That is, the 12 x 12 matrix that was 
developed could have been replaced with two 6X6 matrices. The procedure for devel- 
opii^ this pair of matrices is a fairly obvious modification of the procedure that was 
used for developing the 12 x 12 matrix. Alternatively, the procedure that was used for 
developing the two 9X9 matrices for the mirror with the central obscuration may be 
applied directly to the case of the simple mirror. One has only to put R equal to zero 
in tables DC and XI and recalculate the elements of the [M^ and [M^ matrices. The 
[c J and ^ J matrices, and the [js]] matrices may be used without 

modification. 

The Nonlinear Problem 

Equations (4) and (5) were derived from equation (1) under the assumption that the 
error in the mirror surface was no more than about one-twentieth of a wavelength. 
Results calculated by the exact equation in references 2 and 3 show that as the error in 
the mirror surface increases, the effects on the diffraction pattern quickly acquire a dif- 
ferent character and are not even recognizable as extrapolations of the small-error 
effects. For such cases, experienced opticians can, by observing the pattern as they 
move the knife-edge in the focal plane, locate the error and estimate its magnitude. The 
technique can possibly be adapted for an orbiting telescope, so that it might generally be 
possible to reduce a large error by successive applications of relatively coarse correc- 
tions imtil the error is sufficiently small to be determined exactly by the small-error 
theory. 

One may hope that the problem of reducing large errors would exist only during the 
initial adjustment of the telescope after launch, and that with adequate thermal control of 
the telescope, the errors that develop subsequently will remain within the range that can 
be analyzed by the linear theory. In any case, a reasonably practical method of soMi^ 
the nonlinearized problem, or, at least, of taking into account the second-order terms 
that were deleted from equation (3) would be very desirable. 

Color 

It has been tacitly assumed throughout the present analysis and discussion that the 
light is monochromatic - that is, that the starlight is filtered through a fairly narrow 
band-pass filter. Actually, for a perfect mirror, the diffraction pattern (see fig. 1) is 
independent of wavelength. If the mirror is not perfect, however, the deviation p from 
the ideal intensity distribution of figure 1 is proportional to the error m (see eq. (5)), 
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where m is measured in half-wavelengths. For a given imperfection, then, m will 
be greater for blue light than for red light, and the intensity increment p will accord- 
ingly be greater for blue light than for red light. It follows that use of a blue filter will 
maximize the sensitivity of the knife-edge test. 

On the other hand, if the imperfection in the mirror figure is large, the present 
linearized theory might not be very applicable if blue light is used (see the preceding 
section), but it might remain satisfactory if red or infrared is used (since the imperfec- 
tion will then be smaller in terms of half -wavelengths). Thus, the longer wavelengths 
should be used for analyzing large errors in the mirror figure, but the shorter wave- 
lengths will provide increased sensitivity after the error has been reduced. 

CONCLUDING REMARKS 

It appears from the present work and the work of Gatewood that Linfoot's integral 
equation for the light intensity distribution in the knife-edge pattern can be readily 
inverted to determine errors in the primary optics of a telescope, provided the errors 
are small. This approach to determining the errors seems practical and furthermore 
requires no special apparatus in the main telescope tube or special treatment of the 
mirror (as are required in other suggested approaches). Nevertheless, many details of 
the procedure remain to be studied. Some of these details have been discussed in the 
last section of the present paper. 

Laboratory studies of the method will be needed to identify problem areas and to 
indicate practical experimental approaches to obtaining the basic data. 

Langley Research Center, 

National Aeronautics and Space Administration, 

Hampton, Va., March 29, 1971. 
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APPENDIX 


MATRICES USED IN ANALYSIS 

This appendix contains the four inversion matrices that were used for the calculated 
examples of the present paper. They are applicable to 12-point and 24 -point analyses of 
the simple mirror and to 18-point analyses of the mirror with a central obscuration 
(Cassegrain type). As was mentioned in the main text, the number of significant figures 
shown is doubtless excessive. 

The printouts that are here reproduced are only six columns wide; accordingly, 
each matrix row requires at least two rows of the printout. Thus, in the 12 X 12 matrix, 
the first row of the printout, reading from left to right, contains the first six elements 
of the matrix row; and the second row of the printout, reading from left to right, contains 
the seventh to the twelfth elements. 


38 



12 X 12 INVERSION MATRIX 


FOR SIMPLE MIRROR 


ROW 1 

1.58027O75E-OI 9.A955284CE-02 2.40752492E-02 4. 84739143E-02 2.21414808E-02 3.51095584E-02 

2,01664358E-02 2,71101263E-02 i •799462586-02 2. 03809443E-02 1. 498552 64E-02 1.1 5022 107E-02 

ROW 2 

-4.77232344E-02 1.66533159E-01 1.32233825E-01 I. 27709820E-02 5.84097441E-02 1.43423129E-02 

3.93341090E-02 1.46342319E-02 2 . 826553 04E-02 1 . 3 9090034E-02 1.90056322E-02 1.02244713E-02 

ROW 3 

-J2.* 9^405 54E- 03 -1.04385 188E-01 1.54173409E-01 1 • 61999662E-01 9. 31 196557E-03 6.39332915E-02 

1.13390435E-02 3.98861057E-02 1 •243938836-02 2.59310534E-02 1.25159659E-02 1.24544588E-02 

ROW 4 

-2.03274109E-02 -2.58694095E-03 - 1 ^442 10995E-01 1 • 22214951E-01 1. 85651400E-01 6.49648338E-03 

6.72888099E-02 8.83556405E-03 3^88865613E-02 1. 03592036E-02 2.23813914E-02 9.39169660E-03 

ROW 5 

-6.68016357E-03 -A. 039971746-02 -1.90380095E-03 -1.74884789E-01 7.87 16095 8E-02 1.99043574E-01 

5.18765253E-03 6.628412256-02 7.87504673E-03 3.43213263E-02 1.01439120E-02 1. 36384568E-02 

ROW 6 

^*i6.L9J543£-02 -?• 182921436-03 -5.284222 746-02 -3.03531918E-03 -1.93896871E-01 2.68506993E-02 

2.02459869E-01 3. 58366475E-03 6.20782485E-02 6. 62392926E-03 2.79012469E-02 8.34263162E-03 

ROW 7 

-8.34263161E-03 -2.790124696-02 -6. 623929256-03 -6. 207824866-02 -3. 583664786-03 -2.02459869E-01 

-2.68506993E-02 1.93896871E-01 3.03531 91 8E-03 5.28422274E-02 7. 18292144E-03 1.561985436-02 

ROW 8 

-1.36384568E-02 -1.014391196-02 -3 .432 13263E-02 -7.87504674E-03 -6.62841225E-02 -5. 18765250E-03 

-1.990435746-01 -7.87160958E-02 1.748847896-01 1.90380094E-03 4. 03997 174E-02 6.68016357E-03 

ROW 9 

-9. 39 1696596-03 -2.238139146-02 -1.035920376-02 -3. 888656 1 3E-02 -8.83556403E-03 -6.7288 81006-02 

-6.496483366-03 -1.856514006-01 -1.222149516-01 1.442109956-01 2.586940956-03 2. 032741 C9E-02 

ROW 10 

-1.245445 886-02 -1.2 51596596-02 -2.59 310 5346-02 -1. 243938836-02 -3.988610576-02 -1.133904356-02 

-6.393329156-02 -9.311965566-03 -1.619996626-01 -1.541734096-01 1.04385188E-01 2. 9224 85556-03 

ROW 11 

-1.022608126-02 -1.900681076-02 - 1 .390943486-02 -2. 826509906-02 -1. 4633053 3E-02 -3.933249916-02 

-1.434070306-02 -5.840856566-02 -1.277055076-02 -1.322342566-01 -1.665343376-01 4.772162456-02 

ROW 12 

-1.150221076-02 - 1. 498552646-02 -2.03 8094446-02 -1 . 799462 58E-02 -2. 71 1012636-02 -2.016643586-02 

-3.510955846-02 -2.214148086-02 -4.847391436-02 -2.407524926-02 -9.495528406-02 -1.588270756-01 


00 

CO 
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o 24 X 24 INVERSION [m][f] MATRIX FOR SIMPLE MIRROR 


ROW 1 



6.623402056-02 

1,609706 52E-02 

3.39613915E-02 

1.45984833E-02 

2.53715160E-02 

1.34470894E-02 

2.09345679E-02 

1.25220957E-02 

1.803645146-02 

1.17274715E-02 

1.58738865E-02 




MMIIIIimi IIM 



8.79465651E-03 

9. 53475388 E-03 

7,850034436-03 

7.77355102E-O3 

6.51071174E-03 

4.85198124E-03 





8-77568805 E-02 

8.769272276-03 






2-36080428E-02 

9.456582046-03 



^^rnmFamnatm 

8.98904479E-03 

1. 50 1 63 3 65 E- 0.2 

8.64456724E-03 

1.30731192E-02 

8.20626729E-03 

mSSSwsSm^^^m 

7.6 1871 3 096-03 

9.35606375E-03 

6.726459786-03 

7.14388681E-03 

4.64334420E-03 


ROli_ 3 






-1. 109830 71E-04 

-6. 34986658 E-02 

1.62523356E-01 

1.06912210E-01 

6.62197315E-03 

4.56471843E-02 

7.41519808E-03 

3. 10558516E-O2 

7.79707474E-03 

2.51412410E-O2 

7.96553586E-03 

2.08730335E-02 

8.CC084099E-03 

1.77314273E-02 

7.934746546-03 

1.517398266-02 

7.77167029E-O3 

1.29132411E-02 

7.49045124E-03 

1.07437333E-02 

7.024391966-03 

8.41546916E-03 

6.15416626E-03 

4.97592573E-03 







ROW 4 

-1.07219561E-02 

5.09056917E-04 

-8.90318288E-02 

1.62221895E-01 

1.25745869E-01 

5.02404965E-03 

5.078652f*7E-02 

5.88248521E-03 

3.44215795E-02 

6.40897394E-03 

2.66162330E-02 

6.73563036E-03 

2. 171R4996E-02 

6.92299536E-03 

1,814130146-02 

6.99483119E-03 

1.52407672E-02 

6.94730944E-03 

1.26773006E-02 

6.738859116-03 

1.020397836-02 

6.22681632E-03 

7.48589687E-03 

4.53912092E-03 







ROW 5 






-2.24782096E-03 
4. 177'7l03lE-03 

-2.27914277E-02 

5.52257804E-O2 

1.01936214E-03 

4.986163'43E-03 

-L.11829231E-01 

3.64480075E-02 

1.56281 148E-01 
y55387656E-03 

i.42819297E-0l 

2.76136312E-02 

5.96046962E-03 
7 6_!.^7 83 3 2.27E-oX 

2.2U66702E-02 
L. l94_946'l9E-02 

6.24526321E-03 

6.35832270E-03 

1.81150454E-02 
_9l,028014526-p3 _ ] 

6.42179928E-03 
5 . 85465 149E-03 _ 

1.48614230E-O2 

■_"_5_.d88_04652E-O3 

ROW 6 

-7.a2'32 76 0iE-0y' 
1.58417997E-01 
2.84290058E-02 
1. 4 38 46 316E-a2. 

-1.74'290l03E-03 

3.360"l5‘520E-03 

5.22594541E-03 

5.925559886-03 

-3.16521881E-02 ' 
y.9277'7224E-d2 
2.23631784E-02 
_l.iy)4d525E-02,"' 

7.'lOOOd285E-64 “ 
4.15396303E-O3 ' 
5.58166397E-03 
■ 5.'7469l536E-03 

“*-“171172^97^-01 
3.82395503E-O2 
“ l“.796"26 5 54E-02 * 
" ■7_,85 75’6183E-03‘ 

' l,45d43889E-0l 
4.758277C7E-03 
5.826n373E-03 
4.43419770E-03 

ROW 7 

-3.08196174E-03 
1 .29844603E-01 
4.2'5357B30E-03 
li.524.5J5_l65-03 

-1.52794522E-02 
1.71819574E-01 
2.87352427E-02 
..1-34056765E-02 

-1,044066426-03 
2.87184701E-03 
4.77052372E-03 
7065 5.5 57 E-0 3 

-3.96699966E-02 
6.24757239E-02 
2.21410588E-02 
_ 9*7^2198JE-03 

6.632462106-04 
3.63168524E-03 
5. 19655324E-03 
^ 5^54“9446C9E-p3 ' 

-1.49562970E-01 

3.94261673E-C2 

1.73465295E-02 

5.21353873E-03 

RQW 8 

-6. 72042437E-03 
-1.64813516e-6i 
4.02498795E-02 
U65916654E-02 

-2.9R715125E-03 
i,l09 4“3441E-di 
3.71586193E-03 
5.H533531E-03 

-2.05693919E-02 
l.03i66'392E-O 1 
2.87869939E-02 
U23419350E-02 

-1.105714826-03 

2.320628C9E-03 

4,274489756-03 

5.24970756E-03 

-4.63069365E-02 

6.50496959E-02 

2.17290729E-02 

8.29434092E-03 

3.17077527E-04 

3.07151792E-03 

4.75020908E-03 

4.32000251E-03 

ROW 9 

-3.5331980OE-O3 
1. 701 958fi0F-04 

-1.22336000E-02 

-1.77784924E-01 

-2.376371056-03“ 

a.'93448613F-62 

■ -f.54522i35E-02 " 
1.91932867E-01 

■- 9.T8 y7l 04 Sl- 04“ 
l.’98187855E-03 

■-5i2"2507697E-02 

6.66413320E-02 

2.72508679E-03 

4.03975373E-02 

3. 39762908 E-03 

2.82787122E-02 

4.012893686-03 

2.07978829E-O2 

4^568^8 lJ^8E^g.3_ 

1,52984.8408-02 

_ .5.^7_89998E^0J_ 

,,l._060920_50E-02 

5.21792_004E-0 3 

.y 36073205E-03 

ROW 10 

-6.14581 806E-03 
-5.7Q57'lOOg£-oa 

-3.79665664E-03 
- 1.6996 9763 E-Q4 

-1.593502146-02 

-1.88008500E-01 

-2.41566427E-03 

6.535B7l"76E-02 

-2.931543576-02 

1.9“8259790E-dl 

-1.20572235E-03 

1.55103035E-03 

6.74758683^-02 

1.96695627E-02 

2.30093986E-03 

4.26700859E-03 

4.0116U21E-0? 

1.388255521-02 

3.00372278E-03 

4,707236356-03 

2.748459586-02 

8.83866508E-03 

3.66519248E-03 

4.18861621E-03 

3fl.W_ _ 

82233 96iE-03 

. -1.05936729E-02” 

-3.336997776-03 

'-l.94il2740E-02 

- 2“. Tfi 5”76‘8CE-by" 

■-3 7^78585 31 E-02 

-I.270U419E-03 

-6.1U29U7E-02 

-3.54959688E-04 

-1.95696945E-01 

3.99134796E-02 

2.01746350E-01 

U26887057E-03 

3.56559612S-03 

6.72469534E-02 

1.79476825E-02 

2.05924595E-03 

4..28663696E-03 

3.910854186-02 
1 .17588392E-02 

2. 816286626-03 
4.83884683E-d3 

2.60816049E-02 

5.54484213E-03 

_P0W‘ ' 12 

j -5.79047771E-03 
-3.54341156E-02 

-4.381945RRE-03 

-1.55545923E-03 

-1,338723536-02 

-6,40615.9206-02 

-3.43815311E-03 

-7.00139678E-04 

-2.20310433E-02 

-2.00478559E-01 

-2.46640413E-03 

1.336716486-02 

2.02587096E-Q1 

2.43952^156-OZ' 

9.17328200F-04 

3.33824959E-C3 

6.62040386E-02 

1.60500961E-02 

L. 715253486-03 
4.086118096-03 

3.76531728E-02 

9.56049878E-03 

2.52347991E-03 

4.02894191E-03 
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24 X 24 INVERSION [m] [f] MATRIX FOR SIMPLE MIRROR - Concluded 

13 . 7 1 

- A.0a 894l 91E-03 -9. 560AO8 78E-03 . 9 861 1 8 09E-Q 3 _ - 1 . 60500961E-02 -3 . 3302A959E-O3 -2.^39525156-02 



-2.523479916-03 

-3.765317286-02 

-1.71525348E-03 

-6. 620403866-02 

-9. 173202006-04 

-2.025870966-01 


-1. 336716486-02 

2,004785596-01 

7.00139678E-04 

6.406159206-02 

1.555459236-03 

3.5434U56E-02 


2.466404136-03 

2. 203104336-02 

3. 438153116-03 

1. 338723536-02 

4.381945886-03 

5.790477716-03 

ROW 

14 

-5.54484P136-03 

-4.838846336-03 

-1.175683926-02 

-4,286636966-03 

-l.7947"6a'25*E-02 

-3.565596126-03 


-2.608160496-02 

-2.8162P662E-03 

-3.910854186-02 

-2.059245956-03 

-6.724695346-02 

-1.288870576-03 


-2.017463506-01 

-3.991347966-02 

1,956969456-01 

3.549596886-04 

6.111291176-02 

1.270U419E-03 

-- 

3._27858531E-02 

2.273576806-03 

1.941127406-02 

3.336997776-03 

1.059362296-02 

3,822359656-03 

_ ROW 

15 

-4. 180616216-03 

-8.838665086-03 

-4. 707236356-03 

-1.38825552E-02 

-4. 267008596-03 

-1. 966956276-02 


-3.665192486-03 

-2,748459586-02 

-3-0b372?78E-03 

-4.0116U21E-02 

-2. 300939866-03 

-6.74758683E-02 


-1.55103035E-03 

-1.98259798E-01 

-6.535871766-02 

1.B3008500E-01 

1.699697636-04 

5,70571000E-02 

— 

1.20572235E-03 

2.C3154357E-02 

2.41566427E-03 

1.593502146-02 

3,796656646-03 

6.145818066-03 

ROW 

.1 6 

-5.360732056-03 

-5.217920046-03 

-1.060920506-02 

-5.027899986-03 

-1.529848406-02 

-4.56860178E-O3 


-2.079788296-02 

-4.012893686-03 

-2.82787122E-02 

-3.397629086-03 

-4.03975373E-02 

-2.72508679E-03 


-6.664133206-02 

-1,981878556-03 

-1.919328676-01 

-8,934486136-02 

1.777849246-01 

-1.70195860E-04 


5.22507697E-02 

9. 883780486-04 

2.54522135E-02 

2.376371056-03 

1.223360006-02 

3.533198806-03 

!_R0W 

17 

-4. 320002516-03 

-8.294340926-03 

-5.24970756E-03 

-1. 234193506-02 

-5.115335316-03 

-1.659166546-02 


-4,750209086-05 

-?. 172907296-0? 

-4.274489756-03 

-2.878699396-02 

-3.715861936-03 

-4,024987956-02 


-3.071517926-03 

-6.504969596-02 

-2.320628C9E-03 

-1.831663926-01 

-1.109434416-01 

1, 648135106-01 


-3.17077527E-04 

4.630593656-02 

1.10571_4RZE-03 

2.056939196-02 

2.987151256-03 

6.72042437E-03 

ROW 

18 







-5.21 3S38J36-03 . 

' -5.54944609E-03 

-9.735229826-03 

-5. 706555576-03 

-1.34056765E-02 

■ -5.52457516E-03 


-1.73465295E-02 

-5.196553246-03 

-2.214105886-02 

-4. 770523726-03 

-2.873524276-02 

-4.25357830E-03 


-3.942616736-02 

-3. 631685246-03 

-6.247572396-02 

-2. 871847016-03 

-1. 718195746-01 

-1.29844603E-01 


1.495629706-01 

-6. 632462106-04 

3.96A‘?^966E-02 

1.04406642E-03 

1.527"94522E-02 

3- 08 196 1746-03 








ROW 

19 

-4.43413770B-03 

-7.857561836-03 

-'5.74691536E-03 

-1.116405256-02 

-5/925J5986E~03 " 

■ -1.43846316E-02 


j-^aa 63X3,736 -03 

-1.796.2;655^E-02 

-5.581663976-03 

-2.2'36317846-02 

-5.225945416-03 

■ -2.842900586-02 


-4.758277076-03 

-3.823955836-02 

-4.15396383E-03 

-5.927772246-02 

-'3.360155206-03 

-1.584179976-01 


.-U45041S89X-01. 

. 1.3L724976Er0l __ 

-7.100002856-04 

_ 3.16521881E-Q2 

1.742901036-03 

7.023276016-03 






ROW 

20 






— 

-5,088846526-03 

-1.486142386-02 

-2.761363126-02 

-5.854651496-03 
-6.421799286-03 
-5. 553R7656E-03 

-9.028014526-03 

-1.811504546-02 

-3.644600756-02 

-6. 358322706-03 
-6. 245263216-03 
-4. 986163436-03 

- 1 . f9^^M 9 E-02~ 
-2. 21 166 7"C 26-02 
-5.522578846-02 

'-6'; 4 78 33 2 27 6-0 3 
■"-5.96046962E-03 
"'-4.177710316-03 



-1.42B192_97E-bl 

-1. 56281148 E-01 

l.ll8a9231E-01 

-1.019362146-03 

2.279142776-02 

■■ 2.247820966-03 


ROW 

21 







-4.539120926-03 

-7.485896876-03 

-6.226816326-03 

-1.020397836-02 

-6. 738859116-03 

-1. 267730066-02 


-6.947309446-03 

-1. 524076726-02 

-6. 994831196-03 

-1.814138146-02 

-6.922995366-03 

-2.171849966-02 


-6. 735630366-03 

-2.661623306-02 

-6. 408973946-03 

-3.442157956-02 

-5.882485216-03 

-5.078662676-02 

— - 

-^024049656-03 . 

-1.25745R69E-01 

-1.622218956-01 

0,903182886-02 

-5.090569176-04 

1.C7219561E-02 

_RO^ 

22 






. 

-4. 97 59 2 5 736 -0'3 
-1.Z91324UE-02. 

-6.154166266-03 
-7. 771678296-03 

-8.41546916E-03 

-1.517398266-02 

-77024391966-03" 
-7. 934746 546 -'03 

-l70743 73”33E-02 
-1.773142736-02 

-7.'49045Z24E-03 
-8. 000840996-03 

— 

-2. 087303356-02 
-4.56471843E-02 

-7.96 5 53'5e6E-0 3 
-6,621973156-03 

-2.51412418E-02 

-1.66"9i2210E-bl 

-'7.797074 746-03 
-1. 625233566-01 

-3. 18558516E-02 
6. 349866586-02 

-7.41519B0'8E-O3 

l.l89830-/l6-04 


JiOy_ 23 

-4.6433_4420F-03 
'-8.206267296-03 
-9.*26b34216F-03 
1 1-9.A2 7 54998V-0 3 


-7. 14388681 E-C3 
-1.30731192E-02 
-1.994216836-0? 
-4.061314496-02 


-6W7645978E-03 

-8.64456724e-0'3 

-9.45658284E-03 

-8.76927'227£-0~3 


-9.356063 756-03 _ '-7^618713096-03 ' '-I. 1238 13956-02 

-1. 501633656-02 -8.989b44'796 -03 _ -^722721076-02 

-2.360a0^2JE-02_ -9^^'59_7 78^E:i03' -2. 92J 126 656-02 

-8.7 75 688 056-02 " -1.5 38689396-01" 2".90306796E-d2 


ROW 24 

-4. 85 19'8 1216-03 
■ ' -l.'f0 56'3932E-02 
- 1.587388656-02 
-2.5~3 "7151 60 6-02 


-6.51071174E-03 
-9.581981626-03 
-1. 172747156-02 
-1 .45 984 8 33 E-^2' 


-7. 773551826-03 
-1.253885746-02 
.603645 146-02 
r h 39156-02 


-7. 850034436-03 
-U03006239E-02 
-l._2522095_7E-02 
6Q~970652 E-02 


-9^534753886-03 -l.J9^656 5 1 E-03 

-1.4lO'28747E-p2“" -1.100195636-02 

-2. 093456796-02 -1.344700946-02 

- 1 . 6^40 2 ^ 56-02 " ' - 1 . 36227016 £-0 1 
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9 X 9 INVERSION MATRIX [Ms][cJ[fJ FOR ANTISYMMETRIC AL p, SYMMETRICAL m. 
MIRROR WITH CENTRAL OBSCURATION, R = 5/41. 

ROW 1 



l.57ail803E-0l 

1.20362791E-01 

9. 338466126-02 
3.977403I1E-02 

4.05208065E-02 

1.06904775E+00 

6.07500224E-02 

5 . 928799026-02 

6.094868356-02 

ROW 

2 

■2.956476^96-02 

1.462712336-01 

1.734038206-01 
1. 74876098 E-02 

1.34646196E-01 

1.09126029E4-00 

2.59563586E-02 

9.628417616-02 

3.56134190E-0? 

ROW 

3 

6.916538536-03 

1.069911236-01 

-7.25948690E-02 

5.77587354E-C2 

1.81527792E-01 

l.05415448Ef00 

1.59221076E-01 

4.532387596-02 

9.05423483E-02 

ROW 

4 

-6.241539696-C3 
1. 76864790E-01 

8.42345343E-03 

4.39612621E-03 

-9.58374232E-C2 

1.10757345E+00 

1. 505103416-01 

2.105961426-01 

2.48104787E-02 

ROW 

5 

3. 681792326-03 
9. 489361576-02 

-1.38863744E-02 

8.92280726E-02 

1.64895529E-02 

1.035l4744EfOO 

-1.272909946-01 

I.55178876E-01 

2.107652O2E-OI 

ROW 

6 

-2.149759956-03 

3.14541900E-01 

3.96771461E-03 

-1.37556354E-02 

-1.87096935E-02 

1.13665916E4-00 

8.51756562E-03 

-1.24254270E-01 

9.59834435E-02 

ROW 

7 

2.08315033E-03 
1. 091674616-01 

-7. 99628756 E-03 
2.27672566E-01 

1.08483925E-02 

9,946789i4E-01 

-2.979370086-02 

2 .64467798E-02 

-1.54340971E-01 

ROW 

6 

-2.42245572E-04 

-8.0U95643E-02 

3.36424653E-04 

-1.04C09902E-01 

-2.06092978E-03 

1.23969337E4-00 

2.606040C3E-04 

-8.75U5131E-03 

-2.37526055E-03 

ROM 

9 

7.15900144E-04 
3. 84866309E-02 

-2.3C826048E-03 

-1.36305501E-01 

4.13569866E-03 

3.61485147E-01 

-7.219092226-03 

1 .11625908E-02 

-2.15804235E-02 
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9X9 INVERSION MATRIX |mJ [0^] |f J FOR SYMMETRICAL p, ANTISYMMETRICAL m. 
MIRROR WITH CENTRAL OBSCURATION, R = 5/41. 


ROW 

1 







1.57811803E-01 

5.10618191&-02 

9.054721566-02 

1.033359926-02 

3.686460346-02 

9.352955866-02 

4.995359366-02 

4.208307876-02 

3.50922653E-02 

ROW 

2 







«3,04912253E-02 

6.39975582E-02 

1.734038206-01 

4.685790916-03 

1.26335680E-01 

9.846470456-02 

2.201225366-02 

7. 048490896-02 

2.114760516-02 

HOW 

3 







7.602515516-03 

4.98907800E-02 

-7.737024366-02 

1.649446356-02 

1.815277926-01 

1.013735276-01 

1 .439094666-01 

3.536196036-02 

5.730167866-02 

HOW 

4 







'7.59051850E-03 

9.124837146-02 

9.932748476-03 

1.388998786-03 

-1.06034287E-01 
1 . 173430646-01 

1. 605103416-01 

1 .817903796-01 

1.73724842E-02 

ROW 

5 







5.18702703E-03 
5. 671530246-02 

-2.579926 ?2E-02 
3.265974136-02 

2.113487046-02 

1.275889966-01 

-1.474610C6E-01 

1. 551788766-01 

1.774534886-01 

ROW 

6 







■3. 73372986E-03 
2.31758389E-01 

6. 681791256-03 
-6.207054896-03 

-2.956317566-02 
1. 727159516-01 

1.216434446-02 

-1.531810436-01 

9.598344356-02 

ROW 

7 







4. 91C39670E-03 
1.091674616-01 

-1.827611676-02 

1.394307866-01 

2.326445286-02 

2.051306996-01 

-5.774850056-02 

4.424961976-02 

-2.094711766-01 

ROW 

8 






- 

•9.32403389E-04 

1.30824959E-01 

1.255553906-03 

-1.040C99Q2E-01 

-7.215766856-03 

4.174592176-01 

8.248013676-04 

-2 .390859 006-02 

-5.26388421E-03 

ROW 

9 







8.182776096-03 

1.86621702E-01 

-2.558188766-02 

-4.047749306-01 

4.300595456-02 

3.614851476-01 

-6.785019526-02 

9.056366656-02 

-1.420218796-01 


CO 
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Figure 5 .- Bump at x = 0.6 on simple mirror. Circled points show results calculated 
by 12 -point analysis and by 2^^-point analysis. Ticks on upper curve show input 
p values for 12-point analysis. 
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